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1. Introduction 

The renormalization group analysis of two-nucleon interaction by Birse, McGovern, 
and Richardson [1] reveals that there are two fixed points, where no scales appear and a 
perturbative expansion is possible in the scale free limit. Around one fixed point, so called 
the trivial fixed point at which all the interactions vanish, a usual perturbation theory such 
as QED, chiral perturbation theory (ChPT) for meson sector and one nucleon sector [2, 3] 
can be constructed. Around the other fixed point, so called the non-trivial fixed point, 
the inverse of a scattering length or a binding energy vanishes for two-nucleon systems, 
and an effective range expansion (ERE) [4, 5] has been known for a long time in this limit. 

After S. Weinberg had suggested an application of ChPT to nuclear physics, especially 
to nuclear forces two decades ago [6, 7], a lot of works have been done so far. (For reviews, 
see, e.g., Refs. [8, 9, 10] and references therein.-^ ) Because pions appear as (massless) 
Goldstone bosons due to spontaneous breaking of chiral symmetry of QCD (finite pion 
mass comes out due to explicit chiral symmetry breaking terms, i.e., quark masses), the 
pions are expected to play a dominant role in the long-range interactions. However, 
a difficulty of systematic application of ChPT to the two-nucleon system stems from 
the appearance of small scales such as large ^-wave scattering length and small binding 
momentum of the deuteron, which are much smaller than the pion mass, ~ 140 MeV. 
A fine-tuning is required in reproducing those small scales in the S-wave two-nucleon 
systems, which implies that the low energy constants can be important as well [11, 12]. 

In this work, we revisit the 5*- wave nucleon-nucleon scattering with perturbative pions 
proposed by Kaplan, Savage, and Wise (KSW) [13, 14]. We adopt the limit of the non- 
trivial fixed point as the "fine-tuned" ground state of the two nucleon system, utilize di- 
baryon fields [15, 16, 17, 18], which have the same quantum numbers of two-nucleon ^5*0 
and ^5*1 states and the infinite scattering length or the zero binding energy, as a zeroth 
order two-nucleon core, and incorporate pion clouds pcrturbatively (one-pion-cxchange 
contribution in this work as the first chiral correction"^) around the non-trivial fixed point 
being represented by the di-baryon fields. 

This may not be so unrealistic choice because it was conjectured by Braaten and 
Hammer that "Small increases in the up and down quark masses of QCD would tune 
the theory to the critical renormalization group trajectory for an infrared limit cycle in 
the three-nucleon system. At critical values of the quark masses, the deuteron binding 
energy goes to zero and the triton has infinitely many excited states with an accumulation 
point at the three-nucleon threshold" [19] . Subsequently, a study with an updated chiral 
effective theory reported an estimated critical pion mass, which makes binding energies 
of the deuteron and its isospin partner zero, tuned to be m^"* ^ 198 MeV [20]. Thus the 
physical pion mass is closer to the critical pion mass than the mass in the chiral limit. 

In addition, as shown in Ref. [21], because the expansion around the non-trivial fixed 
point is equivalent to ERE and takes a form of series of the inverse of potential (or the 

^ There have been intensive discussions concerning renormaHzation and counting rules, with a finite or 
infinite ciitofi', in solving the Lippmann-Schwingcr equation with nuclcon-nuckion interactions obtained 
in chiral effective theory. For more details, see a recent review, e.g., Ref. [10], and references therein. 

^See footnote 6. 



2 



on-shell K matrix), when one renormalizes coefficients of the ordinary contact inter- 
actions by using the effective range parameters, one has comphcated nonhnear relations 
of the coefficients. On the other hand, because our expansion scheme matches well to 
ERE wc obtain a linear relation among the coefficients if we include the di-baryon fields 
in the formalism up to the order we considered in this work. Thus the formalism with 
di-baryon fields provides us a way simple and transparent to account for the role of one- 
pion exchange perturbatively in studying the chiral corrections in nuclear systems at low 
energies. 

We calculate the 5*- wave phase shifts of the A^A^ scattering for ^5*0 and ^5*1 channels up 
to one-pion exchange contributions in the framework mentioned above, ^ where D-wave 
mixture in the spin triplet channel is neglected for simplicity.® We also calculate higher 
order terms of the effective range corrections, V2, f 3, and f 4, whose definitions will be given 
below. Those terms have been calculated by Cohen and Hansen [24, 25] with a pionful 
effective field theory (EFT). In this work, wc obtain new corrections to them from the 
effective ranges. After the scattering length and the effective range arc rcnormalized by 
using empirical values, we retain the /x-dependence in the effective range corrections, V2, 
V3, and V4, and study the role of the scale parameter /i in the dimensional regularization 
(DR) and power divergence subtraction (PDS) scheme. 

Usually, in a perturbative calculation, the same couplings appearing in higher orders 
as those being rcnormalized in lower orders by physical quantities, the scattering length 
and the effective range in the present case, are replaced by the rcnormalized couplings, 
and thus explicit /^-dependence can be removed from an amplitude. However, the higher 
order contributions in the S-wave nucleon-nucleon scattering amplitude converge slowly 
and would bring a significant modification to the amplitude, thus, in this work we rather 
retain the /i-dcpcndence in the higher orders. In other words, the rcnormalized amplitude 
is perturbatively scale-independent, and adjusting the scale parameter /j, can be regarded 
as a probe of "optimal higher order contributions". In general, terms corresponding to 
the scattering length, the effective range, and the higher order corrections, V2, Vs, and V4, 
in the calculation are functions of the scale parameter fi. After renormalization of the 
scattering length and the effective range, we find that the higher order corrections, V2, 
V3 and V4, are significantly sensitive to the value of /x, and the sensitivity has a strong 
correlation to the range where the theory has predictive power. 

This work is organized as the following. In section 2, we review pionless effective theory 
with di-baryon fields and briefiy discuss relations between the formalism and the non- 
trivial fixed point. In section 3, we calculate the nucleon-nucleon scattering amplitudes 

'"'Wc basically employ the standard counting rules suggested by KSW [13, 14]. However, as to be 
mentioned, we do not calculate all of the contributions up to next-to-leading order (NLO) in the KSW 
counting rules, and expand the one-pion-exchange corrections around the inverse of the LO amplitude, 
A^^. Thus two-pion-exchange contributions which would be important are counted as next-to-next-to 
leading order (NNLO) corrections. 

^ It is known that the one-pion exchange potential has a non-perturbative short range contribution, 
the strong attractive tensor force, in ^5*1-^^1 channel [12], and the failure of the perturbative treatment 
of the one-pion exchange contribution in the higher order is reported in Ref. [22] . However, this difficulty 
can be avoided by introducing a regulator to the singular part of the one-pion exchange potential [23] . 
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for ^5*0 and ^Si channel with the one-pion-exchange contributions. In section 4, the 
amphtudes are renormahzed by using the empirical values of the scattering lengths and 
the effective ranges. In section 5, the higher order effective range corrections, V2, v^, and 
V4, are obtained, and the fi dependences of the effective range parameters and the phase 
shifts are studied. Finally, in section 6, discussion and conclusions are presented. 

2. Pionless theory with di-baryon fields 

In this section, we briefly review a pionless theory with di-baryon fields. A relevant 
effective Lagrangian reads [16, 17] 



jCn + jCs + , 



(1) 



where jCn is the Lagrangian for standard one nucleon sector, and Cg and Ct are those for 
the di-baryon fields in the ''^yS'o and the ^Si states, respectively. The Lagrangian £jv in the 
heavy-baryon formalism reads 



Art 



iv ■ d + 



{v ■ df - d"^ 
2m N 



N. 



(2) 



where is the velocity vector satisfying the condition v'^ — 1, and m^v is the nucleon 
mass. The Lagrangians with the di-baryon fields read 



iv ■ d + 



[v ■ df - 



Am 



N 



Vs 



sl{N^Pi'^°^N) + h.c. 



IV d+ ^ / + At 



4m 



N 



(3) 
(4) 



where and ti are the di-baryon fields for the ^Sq and ^Si states, respectively, ag^t = ~1) 
A5 1 are mass differences between the di-baryon masses and the two nucleon mass, and 

Us^t are coupling constants of di-baryon-nucleon-nucleon (dNN) vertices. P^^^°'' and P/ 
are projection operators for ^5*0 and ^5"! two-nucleon states, respectively, which read 



Vs 



T2TaCr2 , 



,(^Si) 



Vs 



(5) 



where a and r are Pauli matrices operating in the spin and isospin spaces, respectively. 

Diagrams for dressed di-baryon propagator are shown in Fig. 1. ^ As shown in Ref. [13], 
two-micleon bubble diagrams are summed up to the infinite order, and thus the inverse 
of dressed di-baryon propagator in the center of mass frame is obtained as 



Gs^t + + A,,0 + ^(/i + ip) 



. 2 f^N 



mNVlt 



2 —P + IJ^ + ip 



(6) 



Energy of the nucleon in the center of mass frame reads E = p'^/m^, and we have employed 
DR and PDS scheme where n is the renormalization scale for the loop integration. 
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Figure 1: Diagrams for dressed di-baryon propagator: a double line denotes a di-baryon 
field and a line does a nucleon field. A double line with a filled circle denotes a dressed 
di-baryon field. 




Figure 2: Diagram for S -wave NN scattering amplitude: a double line with a filled circle 
is the dressed di-baryon propagator obtained in Fig. 1. 



A diagram for the S'-wave NN scattering is given in Fig. 2. Corresponding scattering 
amplitude is obtained as 

On the other hand the (S'-wave amphtude in terms of the phase shift is given by 

47r 1 

A = , 8 

mjv pcotOo — ip 

where p cot 5q is expanded in terms of the effective range parameters as 

p cot 5q = 1 — rp^ + • • • . (9) 

a is the scattering length, r is the effective range, and the dots denote higher order terms. 
With only di-baryon fields, Eq. (7) gives the corrections up to the quadratic term of p, i.e. 
p^, but if we include the one-pion-exchange contributions, we have terms in higher orders 
of p. We note that the two nucleon bubble diagram plays a role similar to a kinetic term 
in mass renormalization in a propagator, and the scattering amplitude is normalized by 
the term proportional to —ip. Consequently, with the term proportional to due to 
the normalization of the amplitude, we can easily make linear relations to the effective 
range terms being expanded in the denominator of the amplitude. Equating Eq. (7) with 
Eq. (8), we have 

— = -/^ r = 2-^- (10 

a ruNVit rnj^yit 



'^Diagrams are prepared by using a C++ library, " FeynDiagram" [26] . 
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(I) 



(n-a) 



(n-b) 



(HI) 



Figure 3: Diagrams for contributions from one-pion- exchange to the NN scattering am- 
plitude: a dashed line denotes a "potential" pion field. See the caption in Fig. 1 as well. 



In general, Lagrangian for the di-baryons corresponding to the effective range param- 
eters may be written as^ 



n=0 



iv ■ d + 



Am 



N 



da , 



(11) 



where d = s, t. The terms of the Lagrangian only require the Galilean invariance for two 
nucleons, and are coefficients, with the factor from the leading dNN interaction, 
to match with the effective range parameters. 



3. Di-baryons and perturbative pions 

Effective Lagrangian with di-baryon fields and pions may read 



(12) 



where and Ct^^ arc the chiral Lagrangians for pion and one-nucleon sector, respectively. 
Because wc do not consider pion loops, Ls and Lt takes the same form as given in Eqs. (3) 
and (4), respectively.^ 

Effective chiral Lagrangians and L-^j^ are standard ones in the literature and we 
have 



= -/,2Tr(A ■ A) + ^Tr(x+) , 



(13) 



1 



= N'' [iv D - 2igAS ■ A] N + {v Df - D'^ + ■ ■ ■ N, (14) 

With D^ = d^ + where = d^^], = d^^, and x+ = ^h^^ + ^X^^; X is 
the symmetry breaking term generating the pion mass, x o*^ at leading order, and ^ 
is the chiral field, = exp(iT • tt//^). f^^ is the pion decay constant and gA is the axial 
vector coupling constant. Dots denote terms irrelevant to this work. 

With the one-pion-exchange contributions (so called "potential" pion) at 0{Q'^) shown 
in Fig. 3, we may write down the whole amplitude as 



A — + , 



(15) 



^ A slightly different expression of the Lagrangian for the deuteron pole has been obtained by Griefiham- 
mer. See Eq. (2.1) in Ref. [27]. 

^The chiral Lagrangians with di-baryon fields have been obtained by Soto and Tarriis [18, 28]. 
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where is the amphtude from the tree-level dressed di-baryon contributions without 
pions, As^t in Eq. (7), and A^^ is the amplitude from the one-pion-exchange contributions 
in Fig. 3. We may divide it as 



where the amplitudes, ^{ii), and >l(m) are obtained as 



Am 



A 



/ 9 a 

1 (9A 

Oa-ip \2U 
1 



-1 + ^In 1 + 



4p2 

2 r 



— 2/x H ^ arctan I — ) + 2ip 



m 



4p2 



(III) 



9A 



I -/^^ + + In 



+2ip 



-/X 



2p 



arctan 



(-)' 









,2 



4p2 



with 



Od 



Ana 2 



(16) 



(17) 



(18) 



(19) 



(20) 



The subscripts s, t are dropped from the di-baryon parameters, cr. A, and y for simplic- 
ity. The loop diagrams have been calculated by using DR and PDS scheme, as done in 
Refs. [13, 14]. 

Expanding the one-pion-exchange amplitude A^^ around the inverse of the di-baryon 
amplitude l/A,^, which is associated with the non-trivial fixed point, wc have 



p cot 5q — ip + 



47r 1 
mjv A 



ip + 



An 



1 



rriN Ad + A^ 



~ -|- 



An 1 



An A^ 



rriN Ad rriN Aj 



(21) 



We note that because the one-pion-exchange correction is obtained in the form of A^^/A'^ 
in the above expression, the pole structures due to the di-baryon propagator in and 
74(7/7) disappear, whereas A^^ and A^^ terms appear as interactions (O^ and Od) along 
with y4(7) and ^(//) terms, respectively. 

Inserting Eqs. (7) and (17-19) into Eq. (21), we obtain 



p cot So = Od 
+20d 



9a^n 
16nl 

ml 



u -I- arctan 
2p 



|{V + r.^ + mJln(;^)-p^ 

2p 



+ {oi-p') 



4p2 



(22) 



We note that the right hand side of Eq. (9) is real, and Eq. (22) satisfies the unitary 
condition for the phase shift the phase shift Sq should be real below the threshold 
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of pion production. Satisfying this condition is not so trivial, as recently pointed out in 
Ref. [29]. For example, in the scheme proposed by Kaplan, Savage and Wise [13], one 
performs resummation of the leading contact interaction Cq, and the next-leading contact 
interaction C2 is treated perturbatively in calculation of the amplitude. The unitarity 
would be satisfied order by order in the perturbative expansion. This treatment makes 
the phase shift complex in the higher order corrections, and thus one needs an additional 
expansion of tan^o in terms of the phase shift 5o. 

4. Scattering length and effective range: renormalization 

We renormalize our results with two effective range parameters, scattering length and 
effective range in the expression. 



p cot So = -- + ^rp^ + F{p) , 
a Z 



(23) 



where F{p) is a function for effective range parameters in higher orders of p; we can 
approximate F{p) oc when p <^ 1. 
If we define aadi) and rd{iJ>) as 



arf(/i) 



a 



QArnN 



-/x^ + + mMn ( — 



1 9a^n f . . 



-1 



, (24) 



' 1 , gA'rnN 

—r H 

.2 ^167r/2 



.1 + 



then we obtain the higher order contributions as 

1 2 



n -1 



-jj, + m^) 



(25) 



F(p) = 



9a^n 



167r/2 



+2 



+ 



adifi) 2 ^ 



m 



—rriT^ H arctan 

2p 



2p 



11 v 

2 J 



m 



4p2 



(26) 



For the renormalization of the scattering length a and the effective range r, we use the 
empirical values of the parameters as 



1 



-8.30 (MeV) , ro = 2.73 (fm) , 



(27) 



for the ^5*0 channel and 



36.4 (MeV) , n = 1.76 (fm) , 



(28) 



for the ^Si channel [30]. Conventionally, as discussed in the introduction, the terms, 
l/ad{ij) s-nd rd{fJ^), in the function F(p) are replaced by the physical quantities, 1/a and 
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Figure 4: l/ad{n) '"^(a*) fof^ ^^o o,nd ^Si channels are plotted as functions of the scale 
parameter /i. 



r, respectively. In this work, we rather retain the l/a(i(/i) and r^ifJ^) terms in the higher 
orders as a probe of the optimal high energy and higher order contributions, and study 
the /^-dependence in the higher order terms, the effective range corrections, V2, v^, and 
V4, and the S-wave phase shifts in higher momentum regions below. 

The value of the scale parameter /i has been chosen as /i — [13, 14]. With this 
particular choice of the scale parameter /i, we have 

^ — and — , (29) 



for the ^5*0 and ^5"! channel, respectively. This implies that there are no contributions 
to arf(/i) from the one-pion-exchange corrections, and it is determined solely by the di- 
baryon contribution. It was shown, e.g., in Ref. [31] that if the value of a scale parameter 
is chosen around the mass of the lightest degree of freedom that is integrated out, then 
the results for the phase shifts in the ^Sq channel agree very well with the experimental 
data. Because we consider only the one-pion-exchange in this work, we vary the scale 
parameter /i up to ~ 400 MeV. 

We plot l/ad{iJ>) and r(i{n) as functions of n in Fig. 4. We can see that l/ad{n) vanishes 
for the ^So channel ai /i — 52.6 and 154.6 MeV, whereas there is no point where l/ad{n) 
vanishes for the ^5"! channel. 

For physically meaningful values of /i, i.e. /i > m^, the value of l/ad(/i) increases 
monotonically as /i increases. At sufficiently large /i, we can approximate Eq. (24) as 

ll^, (30) 



ad(l^) 2 

and thus l/ad{ij) increases linearly with respect to /i. For large /i, we can approximate 
TdifJ') as 

1 / A,. ,.2 \ 

(31) 
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Figure 5: S-wave phase shifts for ^Sq channel (left panel) and for "^Si channel (right 
panel): Curves of our results with /i — m^, 2m^, Sm-,^ are labeled by the /i values, and 
those labeled by "Reid93" are from a modern phenomenological potential. 



Because ^'s in the parenthesis are divided by and we have in front of the paren- 
thesis, we have a relatively slow change of rd(/i), and it eventually decreases because of 
the quadratic term in the parenthesis. 

Fig. 5 shows the phase shifts in the ^-S'o (left) and ^Si (right) channels obtained from 
Eq. (23). If one requires that a result from EFT should be insensitive to the choice of 
a value of cutoff, corresponding to the scale parameter fi in this work, the applicable 
range of the present calculation for the 5'-wave phase shifts is quite restricted, less than 
p ~ 50 MeV when we choose the scale parameter /j, as /j, = m„, 2mT^, and StUt^, as 
shown in Fig. 5. However since we integrate out the degrees of freedom heavier than the 
one-pion exchange, it may be reasonable to consider the scale parameter in the range 
rriTr < /i < 2m7r. Within this range, ^5*0 channel shows the best results with the fi close 
to 2m^. This behavior is consistent with the one observed in Ref. [31]. For the 
channel the results with reasonable range of /x show good agreement with each other to 
substantially large momenta. 

It is worthwhile to consider the radius of convergence for which the perturbative ex- 
pansion of Eq. (21) is reasonable. Since the scattering length and effective range terms 
are determined with the empirical values, the term F{p) includes the contributions from 
one-pion exchange and depends on /x. We consider the ratio defined by 

^_ \F{p)\ 

-l/a + rp^/2 

criterion for the convergence. 
Fig. 6 shows the ratio R in the ^5*0 channel with = m^, 2771,^ and 3771,^. Dependence on 
the value is significant, but regardless of the value, one can say that the perturbative 
expansion of the pion-exchange contribution is reasonable for p < 150 MeV. Agreement 
of the phase shifts to the phenomenological ones is good within this range. 
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(32) 




Figure 6: Ratio R in Eq. (32) as a Junction of momentum p with various values of ji. 

5. Higher order effective range corrections 

After renormalization of the scattering length and the effective range, the function 
F{p) in Eq. (26) has no unknown parameters, and depends only on the scale parameter 
Expanding pcot^o to higher orders, one has 



p cot 5o = -- + \rp^ + V2P^ + vzp^ + v^p^ + 
a 2 



(33) 



We renormalize the scattering length a and the effective range r to the physical values, 
and obtain the higher order parameters, ^;2, vz, and V4, ... by expanding F{p) in Eq. (26) 
in powers of p as 



F{p) 



V2P'^ + V3P^ + VaP^ + 



(34) 



We obtain 



V2 = 



^^3 



9a^n 



9Af^N 



16 



32 



167r/2 [ 3 al{ii)m^ 5 ad{fi)ml ml 



1 + 



16 



128 



+ 



16 1 



167r/2 [ ai(^)m6 7 ad{iJ,)ml 3 
16rrf(/x) ^ l_r2(/i)| 



1 + 



5 2 ml 

256 1 



4 r^(/i) 
3 m^ 



9a^n 



7 ml 3 mi 



+ 



512 



5 a\{ix)m% 9 ad{fJ^)m 



- 16- 



1 + 



^rf(/^) 
ad(/^) 



(35) 



(36) 



(37) 



We note that the expressions of V2, V3 and V4 above are obtained in terms of l/ad{n) 
Tdifj)- Those expressions have been obtained by Cohen and Hassen without rd{^J) and 
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in terms of the physical scattering length a [24, 25]. We can reproduce their results by 
choosing // = m^r and neglecting the contribution rain). 

Values of V2, v^, and are shown in Table 1. The values in the first row for each 
channel are obtained from partial wave analysis (PWA). Those in the second row are 
obtained by putting ad{f^) = oF^'^' and rd(/i) = where a^^^' is the physical value of the 
scattering length, a^^^' = oq for the ^5*0 channel and a^^^' = ai for the '^Si channel. These 
results are equal to what Cohen and Hansen have obtained [25]. Those in the third row 
are obtained by putting arf(/i) = a^'^^' and rai/j) = r'^^'^' where r^'^^' is the physical effective 
range, r*''^^' = Tq for the ^5*0 channel and r^^^' = ri for the channel. We find that the 
corrections from the effective range are significant, but it is not enough to reproduce the 
PWA values of V2, fs and V4. Those from the fourth to the eighth row are our results 
with various n values, n — mjr, 2m7r, 3m,r, whereas n — 178 and 330 MeV for ^5*0 channel 
and fj, — 167 and 246 MeV for ^Si channel are fitted so that the values of V2 of PWA 
are reproduced (we find several points to reproduce the value of V2 of PWA as to be seen 
below). It has been recently discussed that the scale parameter may be considered as an 
additional short-range "counter term" provided one allows for a fine tuning of the scale 
value [32]. The results with /i — m^r, 2m7r and 3171^^ show a sensitivity of V2, V3 and V4 to 
the jj, value, whereas the results fitted to the V2 value of PWA fairly reproduce all of the 
values of V2, and ^4 with the adjusted values of the scale parameter yU. 

Values of ^2, fa, and ^4 in the ninth and tenth rows for the ^Sq channel in Table 1 are 
obtained by Fleming, Mehen, and Stewart (FMS) [22] up to next-to leading order (NLO) 
and next-to-next-to leading order (NNLO) in the KSW counting scheme, respectively. 
They employed so called global fitting method where they used the phase shift data in 
the ranges, p = 7 ~ 80 MeV for NLO and p = 7 ~ 200 MeV for NNLO, in fitting the 
parameters. In addition, they obtained slightly different effective ranges, r = 2.65 fm for 
NLO and r = 2.63 fm for NNLO from the empirical value r = 2.73 fm. Their NLO results 
basically agree with those obtained by Cohen and Hansen [24, 25], whereas the NNLO 
results are improved but not adequate to completely reproduce the results of PWA. It 
may indicate the slow convergence of the higher order corrections as discussed by the 
authors [22]. It is also noted that adjusting the yU- value which improves the estimated 
values of ^2, vs, and ^4 in this work seems to play a similar role to the higher order 
corrections. 

Now we plot our results of 1)2, vz, and for ^5*0 and ^5"! channels as functions of the 
scale parameter /x in Fig. 7. We find the values of V2, V3, and are quite sensitive to the 
value of /x: the sensitivity increases as the order of the terms increases with the units we 
have chosen, /m^, /m^, /m^, respectively. Thus ^4, is the most sensitive, ^3 is in middle, 
and V2 is less sensitive to /x than others. 

The phase shifts for the 5'-waves are plotted in Fig. 8 with Eq. (23) and the ji adjusted 
to PWA V2 values. Our results (labeled by "This work") are plotted with // = 178 MeV 
for the ^5*0 channel and = 167 McV for the ^Si channel. We compare our results with 
the ones which are obtained by keeping only two effective range parameters, a and r in 

^°The values of PWA are taken from Table I in Ref. [25]. 
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f2(fm^) 


f3(fm^) 


t>4(fm'') 




PWA 


-0.48 


3.8 


-17.0 




^phy.^ r = 


-3.32 


17.8 


-106.8 




(jphy.^ j.phy. 


0.69 


3.99 


-25.9 




/J, — 


-1.74 


11.1 


-68.0 




fj, = 2m,r 


0.17 


0.71 


-4.85 




H = 3m^ 


-2.75 


17.8 


-117 




= 178 MeV 


-0.48* 


4.40 


-26.8 




H = 330 MeV 


-0.48* 


4.37 


-29.0 




FMS(NLO)[22] 


-3.3 


19 


-117 




FMS(NNLO)[22] 


-1.2 


2.9 


-0.7 




PWA 


0.04 


0.67 


-4.0 




aP^y-, r = 


-0.96 


4.57 


-25.5 




(jPhy.^ j.phy. 


0.44 


0.48 


-2.88 




II — m^r 


-0.05 


1.42 


-8.22 




fj, = 2m^ 


-0.25 


2.39 


-16.2 




fj, = 3m,r 


-3.23 


21.2 


-139 




H = 167 MeV 


0.04* 


0.75 


-4.1 




= 246 MeV 


0.04* 


0.68 


-4.7 



Table 1: Values of V2, v^, and for ^Sq and ^Si channel: values in the first row for 
each of the channels are from partial wave analysis (PWA), those in the second row are 
obtained by putting a^iij) — a^^^^' and r^ilj) — Q in our results, those in the third row by 
putting ad{ij) = a^'^^^' and r^ilj) = r^'^^^', and those from the fourth to the eighth row are 
our results with different values: for the fourth to the sixth row fi = rriT^, 2m.^, and 3m^, 
respectively, and those in the seventh and eighth row, fi value is fixed by using the value 
of V2 of PWA (the values with * above are fitted ones); fi = 178 and 330 MeV for ^Sq 
channel and n = 167 and 246 MeV for ^Si channel. Values in the ninth and tenth rows 
for ^Sq channel are results obtained by Fleming, Mehen, Stewart (FMS) [22] up to NLO 
and NNLO in the KSW counting, respectively. 
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Figure 7: Effective range expansion parameters V2, fa, and for ^Sq channel (left panel) 
and ^Si channel (right panel) are plotted as functions of the scale parameter fi. 



the eflFective range (ER) expansion, and those from a modern phenomenological potential 
(Reid93) [30]. Compared to the curve of ER, the results with one-pion exchange are 
improved in both accuracy and the range where the theory is applicable. Comparing the 
results with the fine-tuned /x values to those with randomly chosen /x's in Fig. 5, we again 
see that both range and accuracy are improved. It is notable that the fine-tuned values 
of n are within the range < n < 2m^, which may support the perturbative treatment 
of pion-exchange interactions in the framework of di-baryon pionless EFT. 

Finally we investigate the contributions of V2, fs, and Vi terms to the phase shifts. In 
Fig. 9 we plot the phase shifts in the ^5*0 channel for ER, and V2, Vs and terms added 
successively with = 178 MeV. The results are compared with those of ReidQS. For the 
momenta < p < 30 MeV, the results of f 2, and V4 almost coincide. At around p ~ 50 
MeV, the variation from V2 to is about half of that from ER to V2, and the variation 
from f3 to V4 is much smaller than the previous ones. At p = 60 MeV, the variation from 
ER to V2 is similar to that from V2 to f 3 , and the variation from to w 4 is almost half of 
V2 and V3. From around p ~ 70 MeV, the variation from V2 to V3 exceeds that from ER 
to V2, and from around p ~ 80 MeV the variation from V3 to V4 exceeds that from V2 to 
^3. Therefore we can say that the effective range expansion up to V4 converges reasonably 
up to p ~ 50 MeV, as previously pointed out in Ref . [33] . 

6. Discussion and conclusions 

In this work, we calculate the phase shifts of the neutron-proton scattering for the ^Sq 
and ^5*1 channels at low energies employing the di-baryon effective theory with perturba- 
tive pions. We include the one-pion-exchange diagrams, and expand their contributions 
around the inverse of the tree-level amplitude from the di-baryon fields to implement the 
expansion of the perturbative pions around the non-trivial fixed point. Loop diagrams 
are calculated by using DR and the PDS scheme, and we study our results by changing 
the value of the renormalization scale parameter /j,. 

We renormalize the scattering length a and the effective range r to the physical values. 
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Figure 8: S-wave phase shifts for ^Sq channel (left panel) and for ^Si channel (right 
panel): Curves labeled by "ER" are obtained from effective range expansion with the first 
two coefficients, scattering length and effective range, those by "This work" are our results 
using /i — 178 MeV for ^Sq and /i — 167 MeV for ^Si channel, and those by "Reid93" 
are from a modern phenomenological potential. 




Figure 9: Phase shifts in the ^So channel for ER, ER+V2, ER+V2+V3 and ER+V2+V3+V4. 
Left panel shows the results for < p < 60 MeV and the right for 60 < p < 100 MeV. 
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The scale dependence associated with the one-pion-exchange contributions are encoded 
in the functions ad{iJ>) and rd{n) which we retain in the higher order terms of the effective 
range expansion as the probe of the optimal higher order corrections. These functions 
take part in the physical observable p cot 60 through the term denoted by F{p). We also 
calculate the coefficients in the higher-order effective range corrections, V2, v^, and V4. We 
obtain new corrections from rd(Ai) to the terms obtained by Cohen and Hansen [24, 25]. 
Those terms do not have unknown constants and are solely determined by the terms 
induced by the one-pion-exchange contributions. The /i dependences of V2, v^, and V4 are 
plotted and we found the sensitivity to fi. More sensitive dependence on fi is observed in 
the higher-order term with the units {fm^, fm^, fm'^); is more sensitive than others, 
f3 is medium, and V2 is less sensitive and exhibits a mild n dependence. By adjusting the 
^ value to the known V2 value, however, we find that we can reasonably reproduce the 
and V4 values obtained from partial wave analysis. If ji value is chosen randomly, it is hard 
to say that all the t;2, "^3 and values are well controlled and determined by the terms up 
to the one-pion-exchange contributions. In addition, this feature can be observed as well 
from the 11 dependence of the phase shifts; applicable range of the present formalism up 
to the one-pion-exchange corrections is narrow, less than p ~ 50 MeV with = ttItt, 2mT^ 
and 3m^. However, if we use our best values of // for i>2, then the range of the applicability 
of the theory is considerably widened, up to p ~ m^. 

Now we briefly discuss the work reported by Soto and Tarrus employing the di-baryon 
fields and the perturbative pions up to NNLO [28]. In that work, the authors included one 
higher order corrections than the present work, whereas they employed a slightly different 
counting rule, especially in the leading amplitude, A^, in Eq. (7), and the different param- 
eter fitting method. Because they focus on the momentum region, p ~ m^, contributions 
from the effective range parameters, 1/a and r, in Ad are regarded as small corrections 
compared to the ip term generated from the nucleon bubble diagram, and are expanded 
as perturbative terms. Thus scattering amplitudes obtained in this expansion scheme 
do not reproduce the pole structures at very low energies in the S-wave nucleon-nucleon 
scattering, and they have fitted parameters by using the phase shifts data at p = 53 ~ 
216 MeV (E' = 3 ~ 50 MeV). In the NNLO calculation, they found similar results to those 
obtained in the KSW scheme without the di-baryon fields up to NNLO [22]; the phase 
shifts for the ^5*0 channel are reproduced well up to p ~ 216 MeV {E ~ 50 MeV), whereas 
the calculated phase shifts for the ^Si channel deviate from those obtained in the accurate 
potential models at higher momenta and are good up to p ~ 137 MeV {E ~ 20 MeV). 
This may imply that the nonperturbative part of the one-pion exchange contribution in 
two-pion exchange diagrams is not properly renormalized in the higher momentum region, 
p > 137 MeV, and the perturbative treatment in the KSW scheme for the ^5"! channel, 
even including the di-baryon field, is stalled [12]. This issue is indeed beyond the scope 
of the present our work including up to the one-pion-exchange contribution, and we may 
come back to it in the future. 

Though the scale parameter has been conventionally fixed as = m„, we present 
that the variation of provides us a simple and useful test to investigate the range of the 
applicability of the theory. What we have obtained above may just imply a confirmation 
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of a previous result; the slow convergence of the corrections from the perturbative pions, 
but we have newly and quantitatively investigated the range of applicability of the theory 
as a function of the scale parameter /i. 
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